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Abstract 

We discuss the non-perturbative renormaUzation group flow of Quantum Elec- 
trodynamics (QED) coupled to Quantum Einstein Gravity (QEG) and explore the 
possibilities for defining its continuum limit at a fixed point that would lead to a 
non-trivial, i. e. interacting field theory. We find two fixed points suitable for the 
Asymptotic Safety construction. In the first case, the fine-structure constant a van- 
ishes at the fixed point and its infrared ( "renormalized" ) value is a free parameter 
not determined by the theory itself. In the second case, the fixed point value of a 
is non-zero, and its infrared value is a computable prediction of the theory. 



1 Introduction 



With the advent of perturbative renormahzation theory in the late fourties of the last 
century Quantum Electrodynamics (QED) matured to a physical theory of unprecedented 
predictive power. With only two input parameters, the electron's mass and charge, it 
is able to describe, or "explain" a wealth of experimental data, often with spectacular 
precision. Later on, in the seventies, the electromagnetic, weak and strong interactions 
were united in the broader but conceptually similar framework of what is now known 
as the standard model of elementary particle physics. Undoubtedly the standard model 
provides a very impressive description of all three interactions but it also highlights the 
very limitations of quantum field theory in its familiar form, namely the fact that there 
is always a set of parameters (masses, couplings, mixing angles, etc.) which, as a matter 
of principle, cannot be computed within the theory and must be extracted from the 
experiment. In QED this concerns only the electron's mass and charge, but in the standard 
model there are already more than two dozens of similar input parameters. In fact, 
while in the construction of the standard model non-abelian gauge fields and spontaneous 
symmetry breaking made their appearance as new ingredients, the way the pertinent 
quantum field theory is "defined" , i. e. how the infinite cutoff limit is taken, remained 
essentially the same as in QED. Using K. Wilson's modern picture of renormahzation 
[5 the procedure of perturbative renormahzation can be viewed as taking the infinite 
cutoff limit at a trivial, or Gaussian fixed point of the renormahzation group flow. The 
dimensionality of its critical manifold decides about the number of undetermined input 
parameters. Since the fixed point is Gaussian, the corresponding scaling dimensions are 
essentially the canonical ones, and as a result the generalized couplings that appear as 
coefficients of field monomials with a mass dimension not larger than 4 are essentially the 
ones which cannot be computed. In QED these are the monomials ipip and ipj'^tpAn with 
canonical dimensions 2 and 4, respectively. 
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If one takes the infinite cutoff limit of some quantum field theory at a non-trivial, 
or non-Gaussian fixed point [2] the corresponding count relies on the actual scaling di- 
mensions given by the renormalization group (RG) flow linearized about this particular 
non-Gaussian flxed point (NGFP). Those scaling dimensions will in general not coincide 
with the canonical mass dimension of some field monomial. Nevertheless, as we shall see 
in a moment, the number s of parameters which are undetermined in principle is again 
given by the stability (attractivity) properties of the fixed point. It equals the dimen- 
sionality of its ultraviolet critical hypersurface =5^uv- By definition, is the set of all 
points in coupling constant space which is mapped onto the fixed point by the inversqi 



Computing the set of all scaling dimensions at a NGFP, provided it exists, is a 
difficult task which requires hard (non-perturbative) computations, in contrast to the 
Gaussian case where (leaving marginal cases aside) pure power counting gives the essential 
picture. 

Therefore, in principle it is conceivable that if QED or the entire standard model, 
say, should possess a NGFP appropriate for taking the infinite cutoff limit the resulting 
number of undetermined parameters, s = dim^uv? turns out smaller than at the triv- 
ial fixed point. If so, the quantum theory based upon the NGFP would have a higher 
degree of predictivity than its "Gaussian" counterpart which is provided by the familiar 
perturbative renormalization procedure. A theory defined by means of a NGFP is some- 
times referred to as "non-perturbatively renormalized" . In this paper we shall encounter 
a version of QED which could possess such a higher degree of predictivity. 

Already in the early years of QED a non-trivial UV fixed point was speculated 
about, for a somewhat different reason though [3]. The well-known one loop formula 



^In this paper the orientation of the RG flow is always from the ultraviolet (UV) towards the infrared 
(IR), i. e. in the direction of the natural coarse graining flow. 



RG flow. 




(1.1) 
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suggests that it might be difficult to obtain an interacting theory in the hmit when the 
UV-cutoff A is removed: Keeping in (11. ip the value ca of the bare charge fixed and sending 
A to infinity one finds that the renormalized charge Cren vanishes, so one is left with a 
"trivial" theory. (Conversely, keeping Cren fixed, it is impossible to let A — )■ oo since ca 
diverges at a finite value of A, the Landau pole.) Clearly, if the exact version of eq. (II. ip 
displays a UV fixed point such that ca e* for A oo the prospects for an interacting, 
cutoff-free theory were much better. 

However, to the best of our knowledge no such fixed point exists. On the contrary, 
comprehensive lattice simulations [4J and studies using non-perturbative functional RG 
methods ^ lead to the conclusion that QED is very likely to be a trivial theory in 4 
dimensions. 

As we shall argue later on, the situation might be different when QED is coupled to 
quantized gravity. 

Trying to take the infinite cutoff limit at a NGFP, provided there exists any, is an 
option even when the theory under consideration is non-renormalizable in perturbation 
theory. A well known example of a perturbatively non-renormalizable theory which can 
be defined in this way (which is "non-perturbatively renormalizable" ) is the Gross- Neveu 
model in 3 dimensions 

The situation seems to be similar in Einstein gravity in 4 dimensions. While the 
perturbative quantization of general relativity leads to a non-renormalizable theory [6H8] 
it now appears quite likely that there exists a NGFP suitable for taking the infinite cutoff 
limit there, leading to a predictive theory with only a small number of free parameters. 
The idea of defining Quantum Einstein Gravity (QEG) at a non-trivial fixed point is 
due to S. Weinberg [9] who coined the term Asymptotic Safety for this scenario. This 
term highlights the analogy with asymptotic freedom, the key difference being that now 
a non-Gaussian rather than Gaussian fixed point controls the ultraviolet behavior. While 
originally the viability of the Asymptotic Safety idea could be tested near 2 dimensions 
only [9j , the gravitational average action introduced later |10j opened the way for detailed 
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non-perturbative studies of QEG in 4 dimensions [TUHT7] . Besides finding significant 
evidence for the very existence of a NGFP it was also realized that, consistent with general 
expectations p], the dimensionality s = dimt5^uv indeed seems to be a small number and 
probably does not increase beyond a certain number when the total dimensionality of the 
truncated coupling constant space is increased. 

This phenomenon was first observed in a 3 parameter RG flow in d dimensions 
which included the running Newton constant, cosmological constant and the prefactor of 
a third invariant, J d'^x^/gR'^, |T3]. The flow equations depend on the real (not necessarily 
integer) number din a continuous way, and it was found that there exists a NGFP for every 
dimensionality in the interval betweerj^ d = 2 + e and = 4 (at least). Interestingly, its 
attractivity properties depend on d in an essential way: Above a certain critical dimension 
(near 3), in particular in ci = 4, the NGFP is UV attractive in all 3 directions, while 
below the critical dimension one of the three directions is UV repulsive. In particular in 
d = 2 + e it was found that s = 2 ! Thus the RG trajectories of all asymptotically safe 
theories are confined to a 2-dimensional surface imbedded in the 3-dimensional parameter 
space. Hence one of the couplings can be predicted in terms of the other two. 

The situation is similar in c? = 4 where, however, the "stabilization" of s at a small 
value is seen only when a larger parameter (or "theory" ) space is used. Making an ansatz 
of the form J d'^Xy/g f{R) where / is a polynomial in the curvature scalar it was found 
that dim^uv stabilized at s = 3 when the degree of the polynomial was increased [T5] . 
In a calculation with 8 free parameters in / this allowed for 5 predictions in terms of 3 
input parameters. 

In this paper we are going to discuss the issue of non-perturbative renormalizability 
and the possibility of enhanced predictivity for QED coupled to quantum gravity. One 
of the motivations are various recent perturbative calculations of quantum gravity cor- 
rections to the QED or Yang-Mills beta function governing the RG running of the gauge 
coupling [T81 - I27] . The to date most complete perturbative analysis uses a gauge fixing 

^Note, however that no e-expansion is performed here [T5] . 
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independent approach and a regulator which retains power-hke divergences [23]. It leads 
to the following 1-loop RG equation for the running of the electric charge with the energy 
scale E: 



The first term on the RHS of (11.21) is the familiar one from the fermion loops which tends 
to increase e at large energies. The second term, the gravity correction, involves Newton's 
constant G and the cosmo logical constant A. It has a negative sign and tries to drive e 
to smaller values as E increases. In fact, it has been claimed [25] that the electric charge 
vanishes at high energies and may be regarded an asymptotically free coupling therefore. 

In the following we shall reconsider this picture in the light of asymptotically safe 
gravity. We shall demonstrate that if QED coupled to QEG is asymptotically safe, there 
exists a second option for the behavior of the electric charge at high energies: it may 
assume a non-zero fixed point value e* 7^ 0. If this option is realized the asymptotic 
behavior of QED + QEG is governed by a non-Gaussian fixed point whose hypersurface 
t5^uv is likely to have a lower dimension than in the corresponding asymptotically free 
case e* = 0. Within a simple truncation of theory space, we find that for the theory with 
e* 7^ the infrared value of the charge, or the fine-structure constant a = e^/(47r), is a 
computable number which is completely fixed by the electron mass in Planck units. 

The remaining sections of this paper are organized as follows. In Section 2 we 
introduce and motivate the RG equations we are going to study, and in Section 3 we 
show that they possess two distinct non-Gaussian fixed points. In Section 4 we solve 
a simplified version of the RG equations analytically and in Section 5 we discuss the 
Asymptotic Safety scenario related to one of the fixed points where the fine-structure 
constant can be predicted. In Section 6 we supplement the investigation by a numerical 
analysis of the RG flow, and the Conclusions will be given in Section 7. 
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2 The RG equations 

We shall use a projected form of the gravitational average action [TO] to describe 
the non-perturbative RG behavior of QED coupled to QEG in terms of a simple 3- 
dimensional theory space, treating the charge e{k), or equivalently the fine-structure con- 
stant a{k) = e(A;)^/(47r), along with Newton's constant and the cosmological constant as 
running quantities. Combining the results of (TU] for pure gravity in the Einstein-Hilbert 
truncation with the findings of [2Sl[2n] for the gravity corrections to the running of a we 
are led to the following "caricature" of the flow equations: 

dtg = f3g = [2 + 7]Nig,X)]g (2.1a) 
dA = Pxi9,X) (2.1b) 
dta = /3a= (^Ah2{a) - ^$^(0)^^ a (2.1c) 



with the coefficient 



A = ^np. (2.2) 

OTT 



Here we consider for illustrative purposes a variant of quantum electrodynamics with np 

"flavors" of electrons. 

Several comments are in order now. 

(A) The equations are written in terms of the dimensionless running couplings g{k) = 
k"^ G{k), X{k) = A{k)/k'^ and a{k) where k is the IR cutoff built into the average action. 
The dimensionless "RG time" is denoted t = ln(fc/fco). 

(B) The first two equations, fl2.1al) and (12.1b|) . are taken to be those of pure gravity in the 
Einstein-Hilbert truncation. The anomalous dimension of Newton's constant, r]i\f{g,X 
and the beta function for the cosmological constant, /3a(5',A), were found in ref. [lOj l 
Neglecting the backreaction of the matter fields on the renormalization in the gravity 
sector is (at least partially) justified by the investigations in [H] where is was found that 



^For the explicit formulae see eqs. (4.41) and (4.43) in ref. |10| . 
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a Maxwell field and one or a few Dirac fields do not qualitatively alter the RG flow of g 
and A; this calculation had assumed free matter fields though. 

(C) For small g the anomalous dimension r]N can be expanded in a power series in the 
Newton constant according to 

Vn = 5i(A)/ii(A, g) = B,{X) {g + i?2(A)/ + . . .) (2.3) 

with functions Bi and B2 given in [lojl. From several non-perturbative calculations [T7] 
we know the function hi{g, A) rather precisely and we find that -Bi(A) < for all A. Those 
calculations show in particular that the running of g{k) does not change very much if one 
approximates A(A;) ~ 0, i?i(A) ~ i?i(0), -B2(A) ~ 0, whence 

r,N{g,X)^B,{0)g<0. (2.4) 

In terms of the standard threshold functions $p(w)|f used in [10] we have explicitly i?i(0) = 
-l/(37r) [24$i(0)-$K0)]. 

(D) The beta function for a in the third equation, eq. f l2.1cp . involves a pure matter 
contribution, written as Ah2{a)a, and the gravity correction oc g taken from |2E]- The 
former has been computed in perturbation theory, the first two terms being (for np = 1) 



/3a(a)|g=o = Ah2{a) a = a 



(2,5) 



To obtain a qualitative understanding it will be sufficient to employ the 1-loop approxi- 
mation 

/12(a) = a. (2.6) 

Indeed, the lattice and flow equation studies mentioned above indicate that there exists 
no non-trivial continuum limit for QED (without gravity), and this means that /3a (a) 1^=0 
has no zero at any a > 0. Therefore /12(a) = /^^/(^a) starts out as /12(a) = a in the 
perturbative regime a < 1, and for larger a it is still known to be an increasing function: 



^See eqs. (4.40) in ref. [ID]. 
^See eqs. (4.32) in ref. [lU]. 
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/i2(«) > 0. To be able to solve the RG equations analytically we shall set /i2(«) = « for 
all values of a. This is a qualitatively reliable approximation since, as we shall see, at 
most a zero of h^ia) could change the general picture. 

(E) As it stands, /3q, applies only above the threshold due to the mass of the electron at 
k = rrie- At < rrie the fermion loops no longer renormalize a. In the full fledged average 
action formalism this decoupling is described by a certain threshold function. Here a 
simplified description will be sufficient where we set A = if /c < mg. 

(F) The gravity contribution on the RHS of fl2.1cp was derived in ref. [23] within a 
truncation of theory space which included the gauge field action ^^^^ 

besides the Einstein-Hilbert terms. Within the approximation considered there the gravity 
correction to dte is seen to be independent of the interactions within the matter sector, 
if any. Therefore it is the same for QED and the non-abelian Yang-Mills field considered 
in [28] so that we may obtain f l2.1cp by simply replacing there the non-abelian gauge boson 
contribution with the corresponding fermion term (of the opposite sign!). In [28] also 
subleading corrections to dfC involving the cosmological constant were found. They, too, 
within their domain of reliability do not change the qualitative picture and are omitted 
therefore. 

(G) Identifying the scales E and k, the two terms inside the brackets on the RHS of 
the perturbative result fll.2p . in our notation, translate to g + ^gX. Thus, for A = 0, the 
perturbative gravitational correction has the same structure as fl2.1cp from the average 
action. Since A is small in the applications below, subleading corrections such as the term 
^g\ are inessential for the qualitative properties of the fiow. 

We close this section with a word of warning. There is considerable debate in the 
literature about the gravitational corrections to the beta function of gauge couplings, 
on their precise form [TSH25] . as well as their usefulness [SB] and observability [27] in 
scattering experiments. We emphasize that within the Asymptotic Safety program these 
are secondary issues which are not (yet) relevant. Our goal is first of all to construct a 
quantum field theory by devising a way to take the infinite cutoff limit of the corresponding 
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functional integral; we do this by replacing the functional integral computation by the 
task of solving an exact RG equation for the effective average action and trying to 
take the continuum limit at a fixed point of its flow. Only once this is achieved one can 
start to analyze and interpret the resulting theory, and only then questions such as those 
above on scattering experiments can (and should) be asked. 

For the time being we are still in the flrst phase, and so the RG equations used in 
this paper should be seen as a tool towards understanding the flow of F^ and the possible 
continuum limits it might hint at. 

In view of the non- universality of the gravitational corrections [281 129] it is also 
important to stress that a priori all our results hold true only for the very deflnition of 
e{k) used here, namely via the prefactor of the j F^-term in F^. Every comparison with 
different definitions in other settings or schemes would require a separate analysis. 



3 The fixed points 

Let us start the analysis of the system (12. ip by finding its fixed points, i. e. common 
zeros ((?*, A*, a*) of all three beta functions. Obviously there is a trivial or Gaussian fixed 
point GFP at g* = \* = a* = 0. 

Furthermore we know that the subsystem of fiow equations for pure gravity in the 
Einstein-Hilbert truncation, eqs. (I2.1a[ I2.1b|) . admits for a NGFP at {Qq^Xq) 7^ 0. This 
fixed point lifts to a NGFP of the full system located at {g^, Aq, a* = 0). We will denote 
it by NGFPi. This fixed point is trivial from the QED perspective, the electromagnetic 
interaction is "switched off" there, while the gravitational selfinteraction is the same as 
in pure gravity. 

There exists a second NGFP, non-trivial also in the QED sense, if the equation 

h,{a) = -^$}(0) (3.1) 

has a solution for some a = a* 7^ 0. We shall see in a moment that this is indeed the 
case. This fixed point we will call NGFP2. 
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In the following we will be particularly interested in an Asymptotic Safety scenario 
with respect to this NGFP. Near the fixed point the (linearized) flow is governed by the 
stability matrix Bij = du^PuXu*) according to 

dtUi{k) = J2 B^j (ujik) - u*) , (3.2) 
j 

where u = [g, A, a). For the system under consideration the stability matrix is of the form 

^dg^g dxPg \ 

B= dgP, d,^, (3.3) 

\ " / u=u* 

Two of its eigenvalues are therefore identical to the case of pure gravity in the Einstein- 
Hilbert truncation, giving rise to the familiar two UV attractive directions p!T|[T2]. The 
third eigenvalue is given by 



da(3a{u*] 



a )a 



(3.4) 



If a* is positive, which will actually turn out to be the case, the sign of daPaiu*) agrees 
with the sign of h'2{ot*)- 

At this point we take advantage of the information from the lattice and flow equa- 
tions studies trying to find a non-trivial continuum limit of QED without gravity. In 
(D) of Section |2] we saw that their negative results suggest that /i2(«*) > holds true 
even beyond perturbation theory. As a consequence, the third eigenvalue da(3a{u*), cor- 
responding to the a direction in the 3-dimensional g-X-a-theoij space, is UV repul- 
sive. With two UV attractive and one repulsive direction the UV critical hypersurface 

pertaining to NGFP2 is a two-dimensional surface in a 3-dimensional space, i. e. 
S2 =dim^uv(NGFP2) = 2. 

In comparison, let us also analyze the eigenvalues of the stability matrix of the other 
fixed point NGFPi. As Pg and /3a do not depend on a in our approximation the first two 
eigenvalues remain the same as for NGFP2. However, for the third eigenvalue we obtain 
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Ah2ia*) - -ni'^)g* + Ah'2{a*)a 



TT 



"==°--$}(0)/<0, (3.5) 

71 
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such that the third direction turns out to be UV attractive as well. Hence, NGFPi has a 
3- dimensional UV critical hypersurface, i.e. si = dim^uv(NGFPi) = 3. The fact that 
S2 < Si reflects the enhanced predictivity of an Asymptotic Safety scenario with respect 
to NGFP2 compared to NGFPi. 



4 Explicit RG trajectories 

Let us now analyze the flow in a simple analytically tractable approximation. For 
that we expand the functions hi and /i2 to first order in g and a, respectively, 

hi{g) = g + Oig"^) and /12(a) = a + O(a^). (4.1) 

Furthermore, we neglect the running of the cosmological constant and fix A = Aq to a 
constant value. The remaining system of flow equations reads 

dtg=[2 + Bi{Xo)g]g (4.2a) 

dta = [Aa - -$}(0) g] a. (4.2b) 



TT 

In this approximation there clearly exists a NGFP2 with fixed point values 

9* = -j^^ and a* = -^$1(0). (4.3) 

In the following we will express the constant i?i(Ao) in terms of the fixed point value g* 
according to -Bi(Ao) = —2/(7*. 

The approximation allows us to solve fl4.2ap in separation; its solution is given by 

9* 

The constant of integration Gq = limk-^Q g{k)/k'^ can be interpreted as the IR value of the 
running Newton constant. The simple RG trajectory (14.41) for g shares a crucial feature 
with any asymptotically safe trajectory of the exact system for pure gravity, namely that 
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it connects the classical regime g{k) ^ Gok"^ for k <^ mpi = Gq and the fixed point 
regime g{k) ^ g* for k ^ mpi. Note that the Planck mass is defined in terms of the 
constant Gq. 

Due to the simplified form of (14.2bp . the RG equation for a is now an ordinary 
differential equation of Riccati type, which can therefore be solved in closed form without 
the need for a specification of the function g{k). Its general solution reads, with $J = 
$1(0), 

1 1 f& ^ ^ 

— = — exp -$ / 

where ao = a(fco) is the value of the fine-structure constant at a fixed reference scale k^. 
If we now specialize for the function g[k) of eq. fl4.4p we can perform the integrations in 
f l4.5p and we find 



a{k) \g* + Gokf 



^ f. . 9* \ f. . . . 9* 



(4.6) 



a* V GoP ; V ' ' Gok\ 
where 2-^1 (0., b, c; z) denotes the (ordinary) hypergeometric function. 

From eq. (14. 6 p we infer that there exist three kinds of possible UV behavior for a{k). 
They differ by the value of the terms inside the square brackets [■ ■ ■ ] on the RHS of (14. 6p . 
This value is independent of the scale k. As the prefactor of [■ ■ ■ ] diverges proportional to 
/^6$;g*/7r ^\^Q^ _s. -\yg gjifi the limit limfc_j.oo a{k) = for every strictly positive value 
[■ ■ •] > 0. This corresponds to asymptotic freedom of the fine-structure constant, and is 
similar to the behavior found by Robinson and Wilczek [18] and Toms [25] , but here with 
a concomitant running of Newton's constant. The corresponding RG trajectories of the 
full system are asymptotically safe with respect to NGFPi. 

For a negative value [• • ■ ] < there will be a scale /clp at which the two terms on 
the RHS of (14. 6 p cancel, so that a diverges at finite energies, corresponding to a Landau 
type singularity. 
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a 

Figure 1. The RG flow on the g'-a-plane imphed by the simplified equations (j4.2p . It is 
dominated by two non-Gaussian fixed points. Their respective value of dim^S^uv differs by 
one unit. (The arrows point in the direction of decreasing k.) 

A third type of limiting behavior is obtained for the case that the bracket vanishes 
exactly: [■■•] = 0. As we are then left only with the second term of the RHS of (14. 6p . and 
since 2-^i(a, c; 0) = 1, we find Imik^oo Oi{k) = a*, corresponding to an asymptotically 
safe trajectory with a non-zero coupling at the FP. This is precisely the behavior to be 
expected due to the UV repulsive direction of the fixed point. Since for A; — oo the 
trajectory will only flow into NGFP2 for one specific value of ao, this value of ao, and 
hence the whole trajectory a{k), can be predicted under the assumption of Asymptotic 
Safety. 

The situation is illustrated by the g-a-phase portrait in Fig.[Tl Bearing in mind 
that the arrows always point towards the IR, we see that NGFPi is IR repulsive in both 
directions shown, while NGFP2 is IR attractive in one direction. This is consistent with 
our earlier discussion which showed that in the 3-dimensional ^f-A-a-space NGFPi has 
3 and NGFP2 has only 2 IR repulsive (or equivalently, UV attractive) eigendirections. 
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In Fig.m the trajectories inside the triangle GFP-NGFP1-NGFP2 are those cor- 
responding to the case [■ • • ] > above; they are asymptotically safe with respect to 
NGFPi. The NGFP2 — )■ GFP boundary of this triangle is the unique trajectory (head- 
ing towards smaller g and a values) which is asymptotically safe with respect to the second 
non-trivial fixed point, NGFP2. 

The diagram in Fig. [1] corresponds to a massless electron for which A keeps its non- 
zero value at arbitrarily small scales. In reality the a-running due to the fermions stops 
near mg, of course. 



5 Asymptotic Safety construction at NGFP2 

Let us investigate the unique asymptotically safe trajectory emanating from NGFP2 
in more detail. First we note that the condition of a vanishing bracket [■ ■ ■ ] in (14.61) is 
self-consistent in the sense, that the resulting function ao(fco) is of identical form as the 
remaining function a{k): 

Second, let us approximate this function for scales k ^ mpi much below the Planck 
scale. Later on we shall need it at A; = mg, for instance, where mg is the mass of the 
electron. Then the argument „^ 2 = g* ( ) ^ 10^^ is extremely large and this will be 



an excellent approximation. Hence we may safely truncate the general series expansion 
of the hypergeometric function. 



2F,{a,a,c,z)- i-z) ^ 



\n{—z) + 2 ilj{n + 1) — ip{a + n) — iIj{c — a — n) 

after its first term, and approximate the resulting factor 1 + G^k^ / g* ~ 1, such that our 
final result for scales k <^ mpi reads 



9* 3.^1 
$1 . 



a{k) a* vr 



9* \ . ..f^^i.* 



(5.3) 
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Here denotes the Digamma function and 7 is Euler's constant. Using fl4.3p in order to 
reexpress the ratio g* /a* we can write (15.31) also in the following form: 



1 A 



a{k) 2 



(5.4) 



Recall that A = -^np is a completely universal constant, sensitive only to the number of 
(hypothetical) electron species. Hence, for k ^ mpi, we recover the logarithmic running 
a{k)^^ = —Alnk + const familiar from pure QED. 

In the opposite extreme of k comparable to, or larger than the Planck mass the 
gravity corrections set in, stop this logarithmic behavior, and cause the coupling to freeze 
at a finite value a{k — )■ 00) = a*. Obviously, along this RG trajectory no Landau pole 
singularity is encountered! 

Note also that according to eq. (15. 4p we have a{k) oc l/np for every value of k. 
As a consequence, if we consider a toy model with a large number of electron flavors, 
all a-values that appear along the RG trajectory can be made as small as we like, and 
this renders perturbation theory in a increasingly precise. At the fixed point we have for 
instance 

a* = 9$}(0)^. (5.5) 

Up 

In an Asymptotic Safety scenario based upon the fixed point NGFP2, within the 
truncation considered, the infrared value of the fine-structure constant oir = \imk^Qa{k) 
is a computable number. Using eq. (15. 4p to calculate aiR we must remember however 
that as it stands it holds true only for k > rrie- When k drops below the electron mass 
the standard QED contribution to the running of a{k) goes to zero, and the gravity 
corrections are zero there anyhow. Hence approximately, dta{k) = for < < rrie. 
Thus eq. (15.41) leads to the following prediction for oir ^ a{me): 



1 _ A 
am 2 



2 In ( ^"j + ln(/) - 7 - {-^\g* 



(5.6) 



As the fixed point coordinates are an output of the RG equations, the only input parameter 
needed to predict oir in this approximation is the electron mass in Planck units, me/mp\. 
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It is tempting to insert numbers into eq. fl5.6p . With mg = 5.11 ■ 10~^ GeV and 
mpi = 1.22 ■ lO^^GeV one finds mg/mpi = 4.19 ■ 10"^'^, and for tfie optimized cutoff [30] 
we fiave ^\ = 1. Tlie value of g* = —2/Bi{\q) depends on tlie value cfiosen for Aq. For 
Ao = or Ao = A* ~ 0.193, tlie fixed point value of A in tlie Einstein- Hilbert truncation, 
we get g* ~ 1.71 or g* ^ 0.83, respectively. From that we obtain 

1 Ao=0 1 Ao=A* 

^ 10.91 nj. or ^ 10.96 n^. (5.7) 

aiR «iR 

We observe that the result is relatively insensitive to the value of g* and/or ${, but it 
scales linearly with the number of electron species, np- 

Obviously, for np = 1, this estimate differs from the fine-structure constant mea- 
sured in real Nature, a ~ 1/137, by a factor of roughly 13. However, even within the 
limits of our crude approximation (14. ip . a serious comparison with experiment must in- 
clude the renormalization effects due to the other particles besides the electron, all those 
of the standard model, and possibly beyond. Within the "rij? flavor QED" considered 
here we could mimic their effect by appropriately choosing np- It would then follow that 
the observed aiR is consistent with Asymptotic Safety at NGFP2 if ni;' = 13. 

It is reassuring that for this large number the value of the natural expansion pa- 
rameter of QED perturbation theory, (a/vr), is rather small already. At NGFP2, for 
example, one has (a*/7r) ^ 0.38 and (a;*/7r) ^ 0.18, respectively. 

Next let us try the full Standard Model (SM) and its minimal supersymmetric 
extension (MSSM)Jj Applying the above discussion to the weak hypercharge rather than 
the electromagnetic U{1) one again has a one loop flow equation of the type dtai = Aaf, 
this time with A = 41/(207r) for the SM and A = 33/(107r) for the MSSM, respectively 
j. Here ai = 5a / {3 cos"^ 9^) where 9w is the Weinberg angle. It is most convenient to 



^For a related discussion see j32j . 
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compare the prediction of Asymptotic Safety to the experimental value at the Z mass. 
From eq. (15.61) with the new value of A and rrie replaced by Mz we obtain (with Aq = 0): 

af^iMz) ^ 1/25.7 (5.8a) 
afss^(M^) ^ 1/41.3 (5.8b) 

As compared to the experimental value a^^^(M^) ^ 1/59.5 both of these predictions are 
too high, the supersymmetric one less so. Clearly we may not take these numbers too 
seriously. After all, while for the reasons discussed above we believe that the one loop 
form of the matter beta functions is a reliable guide with respect to the general structure 
of the RG flow, its quantitative status is questionable. 

Nevertheless the following observation might be of interest. The predictions (15.81) 
turn out larger than the experimental value since in the SM and MSSM the coefficient A 
is too small. As a consequence, the matter driven renormalizations which reduce ai{k) 
when k is lowered are too weak. If we could take the RG equations seriously at the 
quantitative level the conclusion would be that Asymptotic Safety at NGFP2 is possible 
if there exist more particles with a U{1) charge than those of the SM or MSSM. We find 
it remarkable that not very many more seem to be needed; it is sufficient to increase A 
by a small factor of order unity. 

On the other hand, if ultimately it turns out that the standard model coupled to 
QEG is not asymptotically safe with respect to NGFP2 then its RG trajectory would 
be one of those inside the GFP-NGFP1-NGFP2 triangle in Fig.[Tl In this case it is 
asymptotically safe with respect to the other non-trivial fixed point, NGFPi. As for 
being free from divergences and predictive at all energies this is still not too much of a 
drawback, though. It only means that the U{1) coupling is not a prediction but necessarily 
an experimental input. 
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6 Numerical results 



Returning to QED coupled to QEG we shall now go beyond the analytically tractable 
approximation of the previous section and employ exact numerical solutions g{k), X{k) for 
the pure gravity subsystem of eqs. (12. lap and fl2.1bl) . Thereby the exact form of the func- 
tions r]N = Bi{X)hi{g, X) and P\{g,X) as implied by the Einstein-Hilbert truncation [10] 
are used, and then the two coupled equations for g and A are solved numerically as in [T^ . 
Then, for every given RG trajectory k i— )■ {g{k), X{k)), we calculate the corresponding a{k) 
by inserting g{k) into fl2.1cp and solving this decoupled differential equation numerically, 
too. 

Staying within the one-loop approximation of Pa we thus confirm the existence 
of both non-Gaussian fixed points, NGFPi and NGFP2. In accord with the general 
discussion above, the latter is seen to have two UV attractive and one repulsive direction. 
All RG trajectories heading for — ?■ 00 towards NGFP2 lie in its two dimensional UV 
critical surface ^uv- It is visualized in Fig.[2]by a family of trajectories starting on ^^uv 
close to the FP, which were traced down to lower scales k. 

As the backreaction of the matter on the gravity sector is neglected, the flow in a 



projection onto the gf-A-plane is identical to the one of pure gravity (Fig. 2(a)). We can 
therefore classify the trajectories as in [12], being of type la, Ila, and Ilia, when the IR 
value of the cosmological constant is negative, zero, or positive, respectively. 

As we rotate the coordinate frame (Fig. 2(c)|2(b) ), we see how the critical surface 
is bent in coupling space. Especially we note that the fine-structure constant only gets 
renormalized to small values a ^ 1, if the gf-A-projection of the trajectory is sufficiently 
close to the type Ila trajectory of pure gravity (the "separatrix" [12]). This is because 
only these trajectories give rise to a long classical regime with G,A^ const [SHIM]- They 
spend a tremendous amount of renormalization group time close to the Gaussian fixed 
point of the gravity sector. The classical regime of gravity is needed for the logarithmic 
running of a to be of effect. 
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As a concrete example of a trajectory with a long classical regime we consider the 
"realistic" RG trajectory discussed in [31] and [31]. In these references a specific g-X- 
trajectory has been identified which matches the observed values of G and A. It is of 
type Ilia and can be characterized by its turning point (the point of smallest A) whose 
coordinates are 



The turning point is passed at the scale kx ~ 10~^°mpi. To make the numerical solution 
of the RG equations feasible we transform the equations to double logarithmic variables 
using T{k) = \n{k/kT) = \n{k/mp{) + 301n(10) as the RG time variable. The transition 
at the Planck scale between the classical and fixed point scaling regime therefore takes 
place at about r(/c = mp\) = 301n(10) ~ 69. 

Having fixed the (7-A-trajectory to be the "realistic" one, there is a unique asymp- 
totically safe trajectory relative to NGFP2 in the three dimensional coupling space. The 
corresponding a{kT) can now be found by a shooting method: If we start slightly above 

= =^uv(NGFP2) and evolve towards the UV the coupling a{k) will head to infinity 
at a finite scale, while starting below ^uv "will result an asymptotically free trajectory: 
a{k —7- 00) = 0. This trajectory is asymptotically safe with respect to NGFPi, and in 
Fig.[T] it corresponds to one of those inside the triangle GFP-NGFP1-NGFP2. The 
closer we get to ^uv the more the trajectory gets "squeezed" into the corner of the tri- 
angle at NGFP2, ultimately leading to two separate pieces, the NGFPi — )■ NGFP2 
and the NGFP2 — )■ GFP branch, respectively. Due to unavoidable numerical errors any 
starting value q;(/ct) will eventually opt for one of the two cases, but if we fine tune it 
to happen at sufficiently large r, we will end up with a good estimate for the trajectory 
which is asymptotically safe with repsect to NGFP2, the boundary line NGFP2 — ?■ GFP 
of the triangle in Fig.lH 

The result of this procedure is depicted in Fig.|3]where we seinp = 1. As can be seen, 
there is a rapid transition to the fixed point scaling regime at the Planck scale (r ^ 69), 
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Figure 3. Double logarithmic plot of the running couplings of the "realistic" trajectory. 

above which all three dimensionless couplings remain constant at their fixed point values. 
By fine tuning we were able to choose the initial parameter to l/a(/c'r) = 14.65 which 
ensures that the trajectory stays at the fixed point value for about two orders of magnitude 
in k before it shoots away to infinity. 

The scale of the electron mass corresponds to a r- value r(/c = mg) = ln(4.19 ■ 
10^^'^) + 301n(10) ~ 17.55, which is far above the turning point scale. At this scale 
the asymptotically safe trajectory predicts a value a{me)^^ ~ 10.93 which is in perfect 
agreement with (15. 7p . Hence we can conclude that the running of A, as well as the exact 
functional form of hi{g, A), are of little effect to the IR value of a. 

7 Summary and conclusion 

In this paper we coupled Quantum Electrodynamics to quantized gravity and ex- 
plored the possibility of an asymptotically safe UV limit of the combined system. Using 
a simple truncation of the corresponding effective average action we found evidence in- 
dicating that this is indeed possible. There exist two non-trivial fixed points which lend 
themselves for an Asymptotic Safety construction. Using the first one, NGFPi, the fixed 
point value of the fine-structure constant is zero, and its infrared value ctiR is a free param- 
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eter which is not fixed by the theory itself but has to be taken from experiment. Basing 
the theory on the second non-Gaussian fixed point NGFP2 instead, the fixed point value 
a* is non-zero, and the ( "renormalized" ) low energy value of the fine-structure constant 
Q;iR can be predicted in terms of the electron mass in Planck units. In either case the cou- 
pled theory QED + QEG is well behaved in the ultraviolet, there is no Landau singularity 
in particular, and it is not trivial, i.e., the continuum limit is an (electromagnetically and 
gravitationally) interacting theory. 

The key ingredient in the RG equations considered is the quantum gravity contri- 
bution to Pa which was obtained in ref. [28]. It is proportional to ga = G{k)k'^a, and 
its sign is such that for increasing k it counteracts the growth of a{k) caused by the 
fermions. This can lead to two qualitatively different scenarios for the high energy be- 
havior of QED + QEG. In the first one, which is also seen in perturbation theory [TK|I23]. 
the gravitational effects win over the fermionic ones and a{k) is driven to zero in the 
UV: it has become an asymptotically free coupling. For RG trajectories of this type the 
^-dependence of Newton's constant plays no essential role; the decrease of G{k) becomes 
substantial only after a{k) is almost zero already. Instead, in the second scenario, the 
initial (low energy) value of a{k) is such that it has not yet become very small when the 
weakening of gravity due to the decrease of G{k) sets in. In particular in the asymptotic 
scaling regime it decreases rapidly, G{k) = g*/k^, so that the fermions, still trying to 
increase a for k ^ 00, have a better chance to win over the gravitons now. Along certain 
trajectories they indeed do, but what is more interesting is the possibility of an exact 
compensation of the two trends. This is exactly what happens at the second non-trivial 
fixed point, NGFP2, which is characterized by a non-zero a*. 

Note that this second possibility is closely related to the Asymptotic Safety of (pure) 
gravity. It could not be found in perturbation theory which, while using a similar gravity 
correction to Pa, treats the factor of G it contains as a constant and therefore misses the 
weakening of gravity at high scales. 
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The most remarkable feature of the second fixed point is the reduced dimensionahty 
of its UV critical manifold and the resulting higher degree of predictivity than in per- 
turbation theory. We take this as a first hint indicating that after coupling the standard 
model to asymptotically safe gravity it might perhaps be possible to compute some of its 
as to yet free parameters from first principles. 

Acknowledgments: We are grateful to A. Ashtekar and H. Spiesberger for helpful 
discussions. 
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